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Abstract 

We report on the computation of the critical point of the deconfinement phase 
transition, critical indices and the string tension in the compact three dimensional 
U(l) lattice gauge theory at finite temperatures. The critical indices govern the be- 
havior across the deconfinement phase transition in the pure gauge U(l) model and 
are generally expected to coincide with the critical indices of the two-dimensional 
XY model. We studied numerically the U(l) model for N t = 8 on lattices with 
spatial extension ranging from L = 32 to L = 256. Our determination of the infi- 
nite volume critical point on the lattice with Nt = 8 differs substantially from the 
pseudo-critical coupling at L = 32, found earlier in the literature and implicitly 
assumed as the onset value of the deconfined phase. The critical index v computed 
from the scaling of the pseudo-critical couplings with the extension of the spatial 
lattice agrees well with the XY value v = 1/2. On the other hand, the index r\ 
shows large deviation from the expected universal value. The possible reasons of 
such behavior are discussed in details. 
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1 Introduction 



In this article we continue our exploration of the deconfinement phase transition in the 
three-dimensional (3d) U(l) lattice gauge theory (LGT) started in Ref. pp. The partition 
function of the compact version of this model can be written as 



x£A n=0 

where A is an L 2 x N t lattice, S is the Wilson action 



S[u] = p a ^ cosu(p s ) + p t ^ cosu(p t ) (2) 



pt 



and sums run over all space-like (p s ) and time-like (p t ) plaquettes. The plaquette angles 
uj(p) are defined in the standard way. The anisotropic couplings Pt and ft s are defined in 
Ref. ppj. Since we study the theory at finite temperature, periodic boundary conditions 
in the temporal direction are imposed on the gauge fields. 

Let us recapitulate briefly what is known and/or expected about the critical behavior 
of the 3d U(l) LGT at finite temperature. At zero temperature the theory is confining 
at all values of the bare coupling constant [2], while at finite temperature the theory 
undergoes a deconfinement phase transition. It is well known that the partition function 
of the 3d U(l) LGT in the Villain formulation coincides with that of the 2d XY model 
in the leading order of the high-temperature expansion [3]. When combined with the 
universality conjecture by Svetitsky-Yaffe [3], this result makes one to conclude that the 
deconfinement phase transition belongs to the universality class of the 2d XY model, 
which is known to have Berezinskii-Kosterlitz-Thouless (BKT) phase transition of infinite 
order [3 [6]. It is therefore generally expected the critical behavior of the 3d U(l) LGT to 
coincide with that of the XY model. In particular, one might expect the critical behavior 
of the Polyakov loop correlation function T(R) to be governed by the following expressions 

for P > Pc and 

T(R) x exp [-£/£(*)] , (4) 

for P < Pc, t = P c / P — 1. Here, R ^> 1 is the distance between test charges, T is the tem- 
perature and £ ~ e bt is the correlation length. Such behavior of £ defines the so-called 
essential scaling. The critical indices rj(T) and v are known from the renormalization- 
group (RG) analysis of the XY model: T](T C ) = 1/4 and v = 1/2, where T c is the BKT 
critical point. Therefore, the critical indices r] and v should be the same in the finite- 
temperature £7(1) model if the Svetitsky-Yaffe conjecture holds in this case. 

The renormalization-group calculations of the RG flow, presented in Ref. [JJ , gave sup- 
port to the BKT scaling scenario. However, the critical indices have not been computed. 
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The direct numerical check of these predictions was performed on lattices L 2 x N t with 
L = 16, 32 and N t = 4, 6, 8 in Ref. [TJ. Though the authors of Ref. [7J confirm the expected 
BKT nature of the phase transition, the reported critical index is almost three times that 
predicted for the XY model, rj(T c ) ~ 0.78. More recent numerical simulations of Ref. [8] 
have been mostly concentrated on the study of the properties of the high-temperature 
phase. What is important for us here is the derivation of the critical point in Refs. [7J [5] . 
In these papers it was found that, for the isotropic lattice /3 S — (3 t — (3 with L = 32 and 
N t = 8, the pseudo-critical point is (3 pc = 2.30(2) for Ref. [7J and (3 C « 2.346(2) for Ref. [8]. 
Values of (3 above these values were taken implicitly as belonging to the deconfined phase. 
We shall comment on this derivation later since our result for the infinite volume critical 
coupling differs essentially for this choice of N t . 

In our previous paper pQ we have studied the model on extremely anisotropic lattice 
with 0s = 0. In this limit the model exhibits the deconfinement phase transition which 
gives the possibility to study the critical behavior. We presented simple analytical consid- 
eration which showed that in the limits of both small and large (3t such anisotropic model 
reduces to the 2d XY model with some effective couplings. Then we performed numer- 
ical simulations of the effective spin model for the Polyakov loop which can be exactly 
computed in the limit (3 S = 0. We used lattices with N t = 1,4,8 and with the spatial 
extent L £ [64, 256]. Our main goal was to determine the critical index rj supposing that 
the scaling known from the study of the XY model holds also in our case. The main 
conclusion of our investigation was that the value of the index 77 is well compatible with 
the XY value. We may thus assume that at least in the limit (3 S = the 3d U(l) LGT 
does belong to the universality class of the XY model. 

Encouraged by these findings we have decided to simulate directly the isotropic model 
on the lattice with N t = 8. In this paper we present the results of these simulations for 
a number of different quantities. Our general strategy is essentially the same as in the 
previous paper. Namely, we postulate that the scaling laws of the XY model can be 
used to study the critical behavior of the gauge model. We believe that the information 
gathered so far allows for such an approach to be trustworthy. Nevertheless, in doing so we 
have encountered certain surprises. First of all, the infinite volume critical coupling turned 
out to be essentially higher than the values for the pseudo-critical couplings reported in 
Refs. [5] • As a consequence, the values of (3 used in Ref. [8] to study the deconfinement 
phase lie well inside the confinement phase when the thermodynamic limit is considered. 
Secondly, the index v extracted from the scaling of the pseudo-critical couplings with L 
does agree well with the expected XY value v = 1/2. However, the index rj was found 
to be strikingly different from the XY value, namely rj w 0.50. While the value rj « 0.78 
obtained in Ref. [7] could, in principle be attributed to rather small lattices used, L = 32, 
and to an incorrect location of the critical point, our result is almost insensitive to varying 
the spatial extent if L is large enough. 

This paper is organized as follows. In the next section we describe briefly our numerical 
procedure. The result of simulations are presented in the Section 3. Conclusions and 
discussion are given in the Section 4. 
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Figure V. Scatter plot of the Polyakov loop for = 2.00, 2.30, 2.60 on the 32 2 x 8 lattice. 

2 Numerical set-up 

With the aim of calculating the critical indices and then identifying the universality class 
of the 3d U(l) LGT, we simulated the system on lattices of the type L 2 x N t , with N t = 8 
fixed and L increasing towards the thermodynamic limit. In the adopted Monte Carlo 
algorithm a sweep consisted in a mixture of one Metropolis update and five microcanonical 
steps. Measurements were taken every 10 sweeps in order to reduce the autocorrelation 
and the typical statistics per run was about 100k. The error analysis was performed by 
the jackknife method over bins at different blocking levels. 

The observable used as a probe of the two phases of the finite temperature 3d U(l) 
LGT is the Polyakov loop, defined as 

P(x) = l[Uo(x,t) , (5) 

t 

where Uo(x,t) is the temporal link attached at the spatial point x. The effective theory 
for the Polyakov loop is two-dimensional and possesses global U(l) symmetry. Since 
the global symmetry cannot be broken spontaneously in two dimensions owing to the 
Mermin-Wagner-Coleman theorem the expectation value of the Polyakov loop vanishes in 
the thermodynamic limit. On a finite lattice (52gP(x)) = due to U(l) symmetry (if the 
boundary conditions used preserve the symmetry). This is confirmed by the numerical 
analysis on the periodic lattice: in the confined (small (3) phase the values taken by the 
Polyakov loop in a typical Monte Carlo ensemble scatter around the origin of the complex 
plane forming a uniform cloud, whereas in the deconfined (high f3) phase they distribute 
on a ring, the thermal average being equal to zero in both cases (see Fig. [T] for an example 
of this behavior in the case L = 32, where the transition occurs at (3=2. 346(2), according 
to Ref. [8]). What really feels the transition is then the absolute value of P, which has 
been chosen to be the order parameter in this work. It is worth to stress that this kind 
of dynamics is the same presented by the spin magnetization in the 2d XY model. 
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Figure 2: Polyakov loop susceptibility xl on the lattices L 2 x 8, with L = 48, 64, 128. 



3 Results at (3 S = fit 

At finite volume the transition manifests through a peak in the magnetic susceptibility of 
the Polyakov loop, defined as 



Xl = L 2 ((\P\ 2 )-(\P\) 2 ) , P = -l £P(50 



(6) 



The value of the coupling at which this happens is the pseudo-critical coupling, (3 pc . By 
increasing the spatial volume, the position of the peak moves towards the (nonuniversal) 
infinite volume critical coupling, /3 C . In Fig. [5] the behavior of xl around the transition is 
shown for L = 48, 64, 128. The value of (3 pc for a given L is determined by interpolating 
the values of the susceptibility xl around the peak by a Lorentzian function. In Table [T] 
we summarize the resulting values of (3 pc and the peak values of the susceptibility xl 
for the several volumes considered in this work (we included also the determination for 
L = 32, taken from the first paper in Ref. [8]). 

In order to apply the finite size scaling (FSS) program, the location of the infinite 
volume critical coupling /3 C is needed. The scaling law by which /3 C can be extracted from 
the known values of f3 pc {L) depends on the nature of the transition and, in particular, on 
the behavior of the correlation length. There are, in principle, two hypotheses to be tested: 
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Table 1: j3 pc and peak value of the Polyakov loop susceptibility xl on the lattices L 2 x 8. 



L 




XL,max 


32 


2.346(2), Ref. [8] 




48 


2.4238(67) 


12.93(41) 


64 


2.4719(39) 


20.09(66) 


96 


2.5648(96) 


38.8(1.6) 


128 


2.6526(59) 


60.1(3.5) 


150 


2.68(1) 


92.6(8.0) 


200 


2.7336(69) 


144(12) 


256 


2.7780(40) 


220(20) 



first order and BKT transition. The hypothesis of first order transition is not incompatible 
with data for the peak susceptibility for L > 128. However, the corresponding scaling law 
for the pseudo-critical couplings, 

Ppc = Pc+^ , (7) 

seems to be ruled out by our data (x 2 /d.o.f equal to 5.6 for L > 96, 3.7 for L > 128, 2.1 
for L > 96). 

Assuming the essential scaling of the BKT transition, i.e. £ ~ e bt " , the scaling law 
for (3 pc becomes 

fl* = &+ n ,l m i • (8) 
(InL + B)» 

The index v characterizes the universality class of the system. For example, v — 1/2 
holds for the 2d XY universality class. 

We tried at first a 4-parameter fit of the data for P pc {L) given in Table [1] with the 
law given in Eq. (jSJ). We excluded systematically from the fit the data for f3 pc {L) at the 
lowest spatial volumes, looking for a region of stability of the parameters. Defining L min 
as the smallest value of L for which f3 pc (L) has been considered in the fit, we could not 
find a stable fit for L min < 96. In particular, we found that the x 2 /d.o.f. is ~ 10 for 
L m in = 32, 6 for L min = 48 and pa 1.8 for L min = 64. Moreover, we observed a strong 
dependence of the fit parameters on the starting values used in the MINUIT minimization 
code, although the resulting fitting curve turned out to be in general rather stable. The 
instability of parameters becomes less severe when L min increases and, in particular, for 
L m m = 64 the parameter v becomes compatible with the XY value, v = 0.5, although 
still undergoing large fluctuations under change of the starting conditions of the MINUIT 
minimization procedure. A stable fit could be achieved only for L m i n = 96 and these are 
the resulting parameters: 

#. = 3.06(16), A = -5.3(5.1) £ = -1.4(1.0), v = 0.49(16) (xVd-c-.f. = 1.5). 

We repeated then the fit with the law (jHJ) keeping the parameter v fixed at the XY 
value, v = 1/2, thus reducing to three the number of free parameters in the fit. In 
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this case, the fit instability is highly suppressed with respect to the previous 4-parameter 
analysis and, indeed, already for L min > 48 we can quote stable values of the fit parameters 
(see Table [2]). One can see that an acceptable xVd.o.f. and a stable fit are obtained for 
Lrmn = 64 and L min = 96 and that, for the latter volume, (3 C is consistent with the result 
of the 4-parameter fit. We take therefore (3 C = 3.06(11) as our estimation for the infinite 
volume critical coupling. The determination of f5 c is the first main result of this work. 

While performing this work, we considered also the possibility to extract the index 
v from directly fitting the essential scaling law £ ~ e bt against lattice data for the 
correlation length taken for several /3 values and for several volumes. To be more precise, 
for each considered value of L and for several f3 values across (3 pc (L), we determined the 
correlation length £l as the inverse decay length of the 2-point correlator of the Polyakov 
loops, interpolating the latter as if the exponential fall-off with the distance apply even 
above (3 pc , where, in fact, this correlator decays power-like (see Fig. [3]). At each volume 
it happens that £l defined in this way increases with (3 till (3 w (3 pc and then saturates, 
consistently with the fact that the region of power-like behavior has been reached. It 
occurs, however, that the set of all lattice data for the correlation length £l that, at each 
volume, belong to the region (3 < f3 pc (L), lie approximately on the same curve. This is 
expected to occur more and more accurately as the thermodynamic limit is approached. 
One could then try to fit the lattice data for £^ falling on this curve with the essential 
scaling law and extract v. Unfortunately, the quality of our data did not allow us to have 
a stable fit and we had to reject this method. It cannot be excluded, however, that it will 
be reconsidered in possible future studies of the same kind. 

Table 2: Results of the 3-parameter fit of the values of f3 pc (L) with the law ([8]), v = 1/2 
fixed. 





(3 C 


A 


B 


rVd.o.f. 


32 


5.103(50) 


-1523(86) 


20.03(48) 


13. 


48 


4.65(23) 


-699(239) 


13.8(2.1) 


4.3 


64 


3.44(15) 


-42(24) 


2.4(1.4) 


1.2 


96 


3.06(11) 


-4.7(4.3) 


-1.5(1.1) 


0.76 



Once an estimation for f3 c has been achieved, we can use the FSS analysis, which holds 
just at (3 C , to extract other critical indices. An interesting example is the magnetic critical 
index, rj, which enters the scaling law 

Xl(Pc) ~ L 2 ~ n . (9) 

Actually in this law one should consider logarithmic corrections (see J9j [10] and refer- 
ences therein) and, indeed, recent works on the XY universality class generally include 
them. However, taking these corrections into account for extracting critical indices calls 
for very large lattices even in the XY model; for the theory under consideration to be 
computationally tractable, we have no choice but to neglect logarithmic corrections. 
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Figure 3: Correlation length vs (3 on the lattices L 2 x 8, with L=64, 128, 200, 256. The 
correlation length is determined as if the exponential fall-off with distance of the 2-point 
correlator of the Polyakov loop apply for all (3 values. 



Setting the coupling /3 at the value of our best estimation for (3 C , i.e. (3 = 3.06, we 
determined the susceptibilities Xl(Pc) f° r several volumes (see Table [3] for the results). 
Then, following FSS, we fitted the results with the law Xl(/3 c ) = AL 2 ~ V and got 

A = 0.0171(10), r] = 0.496(15) (x7 d -0- f - = 0.60) . (10) 

This is the second main result of our paper. We stress that this value for 77 is by far 
incompatible with the 2d XY value, i]xy — 0.25. The most extreme consequence of this 
finding is that the deconfinement transition in the 3d U(l) LGT at finite temperature does 
not belong to the same universality class as 2d XY spin model. This would contradict the 
Svetitsky-Yaffe conjecture, raising a problem in the understanding of the deconfinement 
mechanism in gauge theories. We will further comment on this issue in the next section, 
discussing possible ways out. 

In such a situation, it becomes particularly useful to have another determination of the 
index 77, by an independent approach. Following the strategy of our previous paper pQ, 
we define an effective 77 index, through the 2-point correlator of Polyakov loops, according 



S 



Table 3: Values of Xl(Pc = 3.06) on the lattices L 2 x 8. 



L 



Xl(/3 c = 3.06) 



48 

64 

96 

128 

150 

200 

256 



5.732(42) 
8.887(76) 
17.16(95) 
25.37(60) 
31.52(75) 
50.1(2.5) 
65.9(4.6) 



to 



hg[T(R)/T(Ro)] 
hg[R /R] 



(11) 



with Ro chosen equal to 10, as in Ref. p. This quantity is constructed in such a way 
that it exhibits a plateau in R if the correlator obeys the law ([3]), valid in the deconfmed 
phase. 

In Fig. H] we present r] e ^ as a function of the distance R for several (3 values on the 
lattice with L = 200. A drastic change in the behavior of rj e g is observed across the value 
/3* ~ 3. In particular for (3 > (3* a plateau develops at short distances, deviations at large 
R being interpreted as a finite volume effect which becomes stronger with increasing j3 
since £ diverges in the deconfmed phase. The appearance of this plateau is an indication 
that the correlator takes the power behavior expected for a BKT transition. 

The analysis of the behavior of r) e s{R) has been repeated setting at our estimated 
value for (3 C , i.e. (3 = 3.06, and increasing the spatial extent of the lattice. It turns out (see 
Fig. [5]) that a plateau develops at small distances when L increases and that the extension 
of this plateau gets larger with L, consistently with the fact that finite volume effects are 
becoming less important. The plateau value of rj e s can be estimated as rj e s(R = 6) on the 
256 2 x 8 lattice and is equal to 0.4782(25); it agrees with our previous determination of 
the index 77. 

The scenario described by Fig. [5] for (3 = (3 C must be valid for any (3 > (3 C , if the 
system undergoes a BKT transition, since the correlator must obey a power law in the 
whole high-/3 phase. We have found that this is indeed the case by performing an analysis 
similar to that shown in Fig. [5] at several (3 values larger than [3 C (see Fig. [6] for the case 
of (3 = 3.50, which leads to r\ « 0.41)). We observe that, in general, rj((3) < r]((3 c ) for 
[3 > f3 c and had we estimated (3 C from the 3-parameter fit with L min = 64, instead of that 
with L min = 96, i.e. 3.44(15) instead of 3.06(11), the resulting r\ would change by little 
and keep still much larger than the XY value 0.25. 
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R 

Figure 4: r/ eS vs R on the lattice 200 2 x 8 for several (3 values. 

4 Discussions 

In this paper we have studied the critical behavior of the 3d U(l) LGT at finite tem- 
perature. We worked on isotropic lattice with the temporal extension N t = 8. The 
pseudo-critical coupling was determined through the peak in the susceptibility of the 
Polyakov loop. The infinite- volume critical coupling has then been computed assuming 
the scaling behavior of the form (jSj). Our fitting gives the value (3 C = 3.06(11). The de- 
confinement phase is the phase where j5 > (3 C . The detailed study of the deconfined phase 
is clearly beyond the scope of the present paper. Nevertheless, this finding has immediate 
impact on the previous studies of the model. A thorough investigation of the deconfine- 
ment phase was performed in Ref. [S|. However, all (3- values used there are smaller than 
the infinite-volume critical coupling. When the thermodynamic limit is approached the 
critical coupling increases so that the numerical results of Ref. [8] would refer rather to 
the confinement phase of the infinite-volume theory. This is indeed the case as we ex- 
plained in the previous section. One sees from Fig. [3] that the correlation length (inverse 
of the string tension) grows till the pseudo-critical value of /3 is reached. Therefore, the 
string tension is non- vanishing for all values of ft used in Ref. [8]. We conclude the much 
larger /3-values are needed than those used in Ref. [8] to really probe the physics of the 
deconfinement phase in the large volume limit. This however might call for very large 
lattice sizes so the feasibility of such a study is not clear at present. 
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Figure 5: r] eS vs R at (3 C on lattices with several values of L. 



Furthermore, the index v has been extracted from the scaling of the pseudo-critical 
couplings (jHJ). Its value does agree well with the expected XY value v = 1/2. Of course, 
it would be desirable to extract this index directly from the correlation of the Polyakov 
loops along the line described in the previous section. 

One of the main results of our paper is the computation of the index r\ which turns 
out to be 7] ~ 0.496. This value is essentially larger than expected and requires some 
discussion. The easiest explanation would be to state that the spatial lattice size used 
(L G [32 — 256]) is still too small to exhibit the correct scaling behavior, hence the wrong 
values for (3 C and rj follow. However, if one makes a plot of (3 pc (L) vs L, one can see, by 
looking at the trend of data, that it is unlikely that /3 C is much larger than our estimate. 
In fact, our fits with the scaling law (jSJ) show that (3 C decreases when L min increases (see 
Table [2]). Therefore, our result is most likely an overestimation. This implies that the 
true f] is likely even larger than what we found. In any case, even if we use for /3 C rather 
unlikely value (3 C = 3.50, Fig. [6] suggests that r\ ps 0.41, much above the XY value. 

The next objection against our result could be the fact that we have neglected loga- 
rithmic corrections to the scaling law ([9]). It looks for us rather strange that logarithmic 
corrections can lead to decreasing rj almost by two times. We want to mention that in- 
cluding naively such corrections into our fits always results in the increasing of rj values, 
though these values are unstable against the maximal lattice size included into fit. Thus, 
although we cannot rule out this possibility, we do not think that neglecting logarithmic 
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Figure 6: r] eS vs R at (3 = 3.50 on lattices with several values of L. 



corrections results in such a wrong prediction for r\. 

Let us give a simple argument why the index r\ can be different from its XY value. 
Consider the anisotropic lattice. We would like to study the limit of large f3 s . In the 
limit (3 S = oo the spatial plaquettes are frozen to unity. That means, the ground state 
is a state where all spatial fields are pure gauge, i.e. U n (x) = V x V* +e , n — 1,2. 
Perform now a change of variables U (x) — > V x U (x)V* +eQ . Then it is easy to see that 
in the leading order of the large- (3 S expansion the partition function factorizes into the 
product of N t independent 2c? XY models. Let us now look at the correlations of the 
Polyakov loops. Since the Polyakov loop is the product of gauge fields in the temporal 
direction, the correlation function factorizes, too, and becomes a product of independent 
XY correlations, i.e. 

IV (1) 08 a = oo,ft) = [T X y(P t )] Nt ■ (12) 
Hence, for asymptotically large R ^> 1, we get 

r u(1) (f3 s = oc,/3 t >/3^) x 

This leads to a simple relation 

7Kft = oo,/3f) = N tVxY . (14) 



RVXY 



(13) 
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Some conclusions could now be drawn. The critical behavior of the 3d U(l) LGT in 
the limit f3 s — > oo is also governed by the 2d XY model. Nevertheless, the effective index 
r] appears to be N t times of its XY value. Now, for (3 S = we have r)((3 s = 0, /3 t cr ) = t\xy- 
This relation and formula ( JT4"]) allow to conjecture that 

Vxy < r){(3 s ,ff r ) < N tVxY . (15) 

P s = corresponds to the lower limit while f3 s = oo corresponds to the upper limit. In 
general, rj could interpolate between two limits with S . Whether this interpolation is 
monotonic or there exists critical value (3^, such that i](f3 s < (3^ r ,(3^ r ) = t\xy an d r\ 
changes monotonically above P^ 1 ", cannot be answered with data we have and requires 
simulations on the anisotropic lattices. In the paper [11] a renormalization group study of 
3d U(l) model at small (3 S will be presented and computations of the leading correction to 
the large (3 S behavior will be given. The results of our computations support the scenario 
that the index rj depends on the ratio (3 s /(3t- Recently, we have obtained the results of 
simulations for N t = 2,4 performed by A. Bazavov [12J. His results also point in the 
direction of our scenario (see Fig. [7]). In Fig. [7] we plot a possible behavior of r\ supposing 
the monotonic dependence. 

Finally, it is worth mentioning that the factorization in the large /5 S limit does not 
affect the index v. It follows from its definition 03]) that in this limit v = 1/2 as in the 
XY model. We expect therefore that v equals 1/2 for all /3 S and is thus universal. 

In view of our results it might be worth to perform numerical simulations for small but 
nonvanishing j3 s and for larger volumes. The feasibility of a study with larger volumes 
and better accuracy relies strongly on the possibility to improve the simulation code. 
Promising directions could be simulations of the dual formulation of the model (possibly 
with a cluster algorithm) or the use of the Liischer-Weisz algorithm [T3] . The development 
of these directions is in progress. 

Acknowledgment 

O.B. thanks for warm hospitality the Dipartimento di Fisica dell'Universita della Cal- 
abria and the INFN Gruppo Collegato di Cosenza during the work on this paper. Numer- 
ical simulations were performed on the linux PC farm "Majorana" of the INFN-Cosenza 
and on the GRID cluster at the ITP-Kiev. Authors would like to thank A. Bazavov 
for interesting discussions and for providing us with the results of his simulations on the 
lattices with Nt = 2,4 prior to publication. 

References 

[1] O. Borisenko, M. Gravina and A. Papa, J. Stat. Mech. 2008 (2008) P08009 
|arXiv:0806.208i1 [hep-lat]]. 



13 




n l i i i i i i i i i i i i i i i i i i i i 

2 4 6 8 10 12 14 16 18 20 

Figure 7: Conjectured behavior of rj in the /3 s -/3 t plane. Results for N t = 4 are taken from 
Ref. [12]. 

[2] A. Polyakov, Nucl. Phys. B 120 (1977) 429; M. Gopfert, G. Mack, Commun. Math. 
Phys. 81 (1981) 97; 82 (1982) 545. 

[3] N. Parga, Phys. Lett. B 107 (1981) 442. 

[4] B. Svetitsky, L. Yaffe, Nucl. Phys B 210 (1982) 423. 

[5] V. L. Berezinskii, Sov. Phys. JETP 32 (1971) 493. 

[6] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6 (1973) 1181. 

[7] P. D. Coddington, A. J. G. Hey A. A. Middleton and J. S. Townsend, Phys. Lett. B 
175 (1986) 64. 

[8] M . N. Chernodub, E M. Ilgenfritz, A. Schiller, Phys. Re v. D 64 (2001) 054507 
|arXiv:hep-lat/0105021| ; Phys. Rev. Lett. 88 (2002) 23 1601 |arXiv:hep-lat/0112048] ; 
Phys. Rev. D 67 (2003) 034502 |arXiv:hep-lat/0208013| . 

[9] R. Kenna and A. C. Irving, Nucl. Phys. B 485 (1997) 583. 



14 



[10] M. Hasenbusch, J. Phys. A 38 (2005) 5869. 

[11] O. Borisenko, V. Chelnokov, Renormalization group study of 3d U(l) lattice gauge 
theory in the finite temperature limit, in preparation. 

[12] A. Bazavov, private communication. 

[13] M. Luscher and P. Weisz, JHEP 0109 (2001) 010 |arXiv:hep-lat/0108014] . 



15 



